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Abstract 

The width of the neutrino decay into the electron and W boson in a strong external 
magnetic field is obtained from the imaginary part of the neutrino self-energy. This result 
corrects the formulae existing in the literature. The mean free path of an ultra-high energy 
neutrino in a strong magnetic field is calculated. An energy cutoff for neutrinos propagating 
in a strong field is defined. 

1 Introduction 

The first and unique registration of extragalactic neutrinos from the supernova SN1987A ex- 
plosion in the Large Magellanic Cloud, a satellite galaxy of our Milky Way, was undoubtedly 
an exciting achievement of neutrino astrophysics. The solution to the solar- neutrino puzzle 
in an experiment at the heavy-water detector installed at the Sudbury Neutrino Observatory 
was one more important result in this field. This experiment confirmed B. Pontecorvo's key 
idea concerning neutrino oscillations and, along with experiments that studied atmospheric and 
reactor neutrinos, thereby proved the existence of a nonzero neutrino mass and the existence of 
mixing in the lepton sector, see e.g. [HE] and the references cited therein. In this connection, the 
problem of studying the possible effect of an active environment, including a strong magnetic 
field, on the dispersion properties of the neutrinos becomes quite important. 

An analysis of the effect of an external medium on neutrino properties relies on calculating 
the neutrino self-energy operator from which one can extract the neutrino dispersion 

relation, and in part the imaginary part of the neutrino self-energy in medium, denning the 
width of the neutrino decay into the W + boson and a charged lepton, v — > l~W + . Here, we 
consider an electron as a charged lepton, but all the formulas are valid to the muon and r lepton 
as well. 

A literature search reveals that calculations of the neutrino dispersion relation in external 
magnetic fields have a long history [3HlO| . To compare the different results we analyse the 
neutrino self-energy operator that is defined in terms of the invariant amplitude for the 
transition u e — > v e by the relation 

M{v e ^u e ) = - [u e (p) E(p) u e (p)} = -Tr [E(p) p(p)\ , (1) 

where p = (E, p) is the neutrino four-momentum, p{p) = v e {p)v e (p) is the neutrino density 
matrix. On the other hand, the additional energy AE acquired by a neutrino in an external 
magnetic field is defined via the invariant amplitude ([!]) as follows: 
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AE = -^M(u e ^u e ). (2) 

The 5 matrix element for the transition v e — > u e corresponds to the Feynman diagrams 
shown in Fig. [I] where double lines denote exact propagators in the presence of an external 
magnetic field. A detailed description of the calculational techniques for the neutrino self- 
energy operator in external electromagnetic fields can be found e.g. in Ref. [5], see also 
[8llll|. The relevant 5-matrix element can be used to deduce, in a standard way, the invariant 
amplitude (pQ), whence the neutrino self-energy operator takes the form 



7 Q L J2° 0) ^ L + 4p (m e R - m u L) j(*> (p) (m e L - m v R) 
m w 



(3) 



Here, g is the Standard Model electroweak coupling constant; 7 Q are the Dirac matrices; L = 
(1 — 75)/2 and R = (l + 7s)/2 are, respectively, the left- and the right-hand projection operator. 
The integrals introduced in ([3]) have the form 



J m (P) = / (0S( 9 )D<*)(,-p), (4) 

where S(q), G^\q — p) and D^(q-p) are the Fourier transforms of the translation-invariant 
parts of the propagators for the electron, the W~ boson, and the charged scalar boson, 
respectively. We note that the quantity m u in ([3]) is in general the nondiagonal Dirac neutrino 
mass matrix with allowance for mixing in the lepton sector. By this means the flavor non- 
conserving decays u e — > fiW, tW are also possible which are suppressed, however, by the very 
small parameter ~ (m u /mw) 2 < 10~ 22 . 
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Figure 1: Feynman diagrams representing the magnetic-field- induced contribution to the neu- 
trino self-energy operator in the Feynman gauge. Double lines correspond to the exact propa- 
gators for the charged lepton, the W boson, and the nonphysical scalar charged $ boson in an 
external magnetic field. 

The general Lorentz structure of the operator in a magnetic field, defined in Eq. (|3|), 
can be represented in the form [8] 

E(p) = [Al {pi) + &l (pj)\\ + c l (wt)] l 

+ [Ar (jry) + Br (pt)|| + Cr ip<Pl)] R + m u [Ki + i/C 2 (7^7)] • ( 5 ) 

The Lorentz indices of four- vectors and tensors within parentheses are contracted consecutively, 
e.g. {pipy) = p a (p a /3 r y l3 . Further, ip is the dimensionless tensor of the electromagnetic field, 
normalized to the external -B-field, whereas (p is its dual, 

Vap = -Jr 1 <PaP = \ ^al3^ U ■ (6) 
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Finally, in the frame where only an external magnetic field B is present, we take the spatial 3-axis 
to be directed along B. Four- vectors with the indices _L and || belong to the Euclidean {1,2}- 
subspace and the Minkowski {0, 3}-subspace, correspondingly. For example, p± = (0,pi,p2,0) 
and p\\ = (po,0,Q,p3). For any four-vectors P and Q we use the notation 

(PQ)H = (P^Q) = P Q -P 3 Q 3 , 
(PQ)± = (p<p<pQ)=PiQi + ViQ2, 
(PQ) = (PQ)\\-(PQ)x. (7) 

The coefficients Ar, Br, Cr, and K.\^ in ([5]) stem from the Feynman diagram involving 
the scalar <I> boson, while the coefficients Al, Bl, and Cl contain the contributions from both 
diagrams. We note that the coefficients Al, Ar, and IC\ in (|5j) contain an ultraviolet divergence 
which is removed by the vacuum renormalization of the neutrino wave function and mass. 

Using Eqs. (HJ, ([2]) and ([5]), the neutrino additional energy AE in an external magnetic field 
can be written in the form: 

AE = B L ^ [1 - (sv)] + B R i [1 + (sv)] 



m 



[C L - C R + 4/C 2 - (B L - Br) (bv)] (sb t ) + -f (sb 
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+ ^(A l + Ar + 2JC 1 ) , (8) 

where v = p/E is the neutrino velocity vector, s is the unit vector of the doubled neutrino 
spin, b is the unit vector along the magnetic field direction, and h t ^ are its transversal and 
longitudinal components with respect to the neurino momentum, b = b^ + b^ . 

In the previous papers, the neutrino self-energy operator ([3]) was calculated in different 
regions of values of the physical parameters, however, the list of these considered regions appears 
not to be comprehensive. Namely, the investigated limiting cases were the following: 

i) a weak field case (eB <C ml) 

ii) a case of a moderately strong field (rr? e <C eB <C m w ) an( ^ limited region of a neutrino 
transverse momentum with respect to the magnetic field (p± mw) [7]; 

iii) the situation where the neutrino transverse momentum p±_ is rather high, for example, 
P± ^ m w or p± 3> mw, while the magnetic field strength is not too high, eB <C m\, 
which corresponds to the crossed-field approximation 



There is yet another region of values of the physical parameters that requires a dedicated 
analysis. We mean here the case of the high neutrino transverse momentum, when the magnetic 
field strength is also rather high, thus, the crossed-field approximation is not valid. 

This region of parameter values is of importance in connection with problems of the physics 
of magnetars, the pulsars with superstrong surface magnetic fields (B s ~ 10 15 G). In particular, 
the possibility of detecting cosmic neutrinos of ultrahigh energy, ~ 1 PeV or even higher, from 
magnetars is widely discussed (see, for example, [T2TTH] ). It looks reasonable that the process 
of emission of neutrinos having such energies cannot be described adequately without taking 
into account their interaction with a strong magnetic field of a magnetar. 

This talk is based on our recent paper |15] . 



2 Charged-lepton, W- and $-boson propagators in a magnetic 
field 

The Fourier transforms of the translation-invariant parts of the exact propagators in an external 
magnetic field, entering into expressions can be presented in the Fock proper-time formalism 
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in the following form, see e.g. [5j[16]. The lepton propagator is 



S( 9 ) = e-«. { [(„), + m,] \cosf,s - ^ sin J _ Mkl 

Jo cos ps I L ^ cos ps J 



(9) 



where /3 = el? and m e is the electron mass. 

Similarly, the VF-boson propagator can be written as 



G P a(q) = ~ T-e 

J cos ps 



00 rU 



(10) 



- (w) P a cos 2/3s - <^ pCT sin 2/3s 
And finally, for the <l?-boson propagator one obtains 

roo 

D^(q) = / dse" 1 ^, (11 



JO 

where we have chosen the Feynman gauge for the W and $ bosons and have introduced the 
notation (J = e, W) 

n^s^-gD+^ql. (12) 

3 The neutrino decay v — > e~W + in an external electromagnetic 
field 

The probability of the neutrino decay v — > e~W + in an external electromagnetic field is one 
of the most interesting results that can be extracted from the neutrino self-energy operator. 
This probability can be expressed in terms of the imaginary part of the amplitude ([T|) with the 
neutrino self-energy operator ([5]). 

For simplicity, hereafter we neglect the neutrino mass m u , taking the density matrix of the 
left-handed neutrino as p(p) = (jry) L. One obtains: 

w(v -> e~W + ) = ^JmM(u e -»• v e ) = -^ImTr [S(p) (pj) L] = -2^Imfi L . (13) 
E EE 

An analysis of the neutrino decay v — > e~W + in an external field is of interest only at 
ultrahigh neutrino energies. 

In all previous papers the neutrino decay width in an external electromagnetic field was 
calculated in the crossed field approximation, in which case the width is expressed in terms of 
the dynamical field parameter x an d the lepton mass parameter A: 

X = <P FF f /2 , A = 4_. ( 14 ) 

The particular case of a crossed field is in fact more general than it may seem at first glance. 
Really, the situation is possible when the field dynamical parameter x of the relativistic particle 
propagating in a relatively weak electromagnetic field, F < B e (where F means the electric 
and/or magnetic field strength, B e = m^/e ~ 4.41 x 10 13 G is the critical field value), could 
appear rather high. In this case the field in the particle rest frame can exceed essentially the 
critical value and is very close to the crossed field. Even in a magnetic field whose strength is 
much greater than the critical value, the result obtained in a crossed field will correctly describe 
the leading contribution to the probability of a process in a pure magnetic field, provided that 
X S> B/B e . In the frame where the field is pure magnetic one, the dynamical field parameter 
takes the form: 

X = 3— ■ (15) 
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A general expression for the decay width can be written in this case in the form [8] 



\ ni{ v f dz 
o 

x [z + X(l-z)}[2(l + z)(2 + z) + \(l-z)(2-z)}, (16) 
where K 2 /^{U) is the modified Bessel function, with the argument: 

2 [, + A(l-,)]3/3 

Taking in Eq. ()16p the limit A < 1, one obtains the result which can be written in terms 
of the only modified dynamical field parameter 

V A m e mfa 

The range for the £ parameter appears to be rather large, < £ <C l/\/A = m\y/m e . For the 
electron, it is: < £ < 1.6 x 10 5 . If one replaces m e — > m T , then the range is not too wide, 
< f < 45. 

The decay width takes the form 

W( ^ e w ) = ^^^ (, 1 + Tj ex H"Tj • (19) 

The formula (fT9|) should be compared with the results of Refs. [4.6.9J. It should be mentioned 
that the decay width w defined in Refs. [HE] is the same, in the natural system of units, 
than the absorption coefficient a |6j and the damping rate of the neutrino 7 [9j. One can 
see that the absorption coefficient a presented in Eq.(25) of Ref. [6j looks very similar to our 
Eq. (|19p . However, the angular dependence in our formulas is quite different: instead of the 
factor p\/E = Esm 2 6 standing in our Eq. f)19[) . there is the factor p± = Esin9 in Eq.(25) of 
Ref. 0. 

On the other hand, one can see that our result (fl"9j) surely contradicts the Eq. (58) of Ref. [9], 
where an attempt was made of reinvestigation of the process v —> e~W + in the crossed field 
approximation. The difference is the most essential at small values of £, where the result of 
Ref. [9] appears to be strongly underestimated. 

In the earlier paper by Borisov et al. [4] the calculations of the process v — > e~W + width 
were performed in the two limiting cases of the small and large values of the parameter \- I n 
the limit x 1 \ (i.e. £ <C 1) their result can be presented in the form 



V2G F 
w = — — — m e 



TT 



eBsinflexpf-v^ Il^k) , ( 20 ) 

V eB p± J 



and can be reproduced from the general formula (|19|) . 

On the other hand, in the limit x ^ 1 (£ ^ l/V^A) the result of Ref. [3] can be written as 

\/3 G*f 

w = — -= — my/ eB sin 9 , (21) 

V2 7T 



and can be reproduced from our more general formula (|16p , 

A problem of the decay v — > e~W + has a physical meaning only in the fields of the pulsar 
type, where the field strength is of order of the critical value ~ 10 13 G. The above formulas 
for the probability except for Eq. (|2ip are applicable for relatively weak fields only, B <C 10 13 
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G. Taking into account the discovery of magnetars which are the neutron stars with the fields 
~ 10 14 — 10 15 G, it is interesting to calculate the probability of the process v — > e~W + in such 
fields where the crossed-field approximation is inapplicable. 

Thus, we will use the following hierarchy of the physical parameters: p\ 3> rriyy ^> eB 3> 
m 2 . A general expression for the process v — > e~W + probability can be obtained by the 
substitution of Eq. © into Eq. (fT5|) with taking account of Eqs. Q - ([TT]) . After calculations 
which are not difficult but rather cumbersome, the process width can be expressed via the 
function §{rj) depending on the one parameter 7] only: 

V = A ■ (22) 



m 



w 



The process width can be presented in the form 



w(u -> e-W+) = ° F {eB l /2 P± , (23) 



where 



^ 1 f dy (tanhy) 1 " (sinhy) 2 — y tanhy 

i] J y 1 / 2 (sinhy) 2 (y — tanh y) 3 / 2 



o 



y tanhy 
r/(y — tanh y) 



(24) 



We stress that we have obtained this formula neglecting the electron mass as the smallest 
parameter in the hierarchy used. 

The formulas (j23|) . (|24p present our main results and are valid in a wide region of the 
parameter r\ values, < r] <^ m 2 v /{eB). The function $(//) is essentially simplified at large and 
small values of the argument. 

In the limit r\ 3> 1 one obtains: 

1 



<%» 1) - 3^(1-0.3), (25) 

and the error is less than 1 % for r/ > 10. 

The formulas (|23|) . (|25p reproduce the probability (|19|). where the limit £ 3> 1 should be 
taken. 

In the other limit 77 1 one obtains 

$(7?<l)~exp^-i^l-i7?+^ 2 ^ (26) 

and the error is less than 1 % for r] < 0.5. 

The formulas obtained allow to establish an upper limit on the energy spectrum of neutrinos 
propagating in a strong magnetic field. Let us take the typical size R of the region with the 
strong magnetic field as R ~ 10 km. If the neutrino mean free path A = 1/w is much less than 
the field size, A <C R, all the neutrinos are decaying inside such the field. For A = 1 km <C R, 
we can find the cutoff energies E c for the neutrino spectrum, depending on the magnetic field 
strength. This dependence is calculated numerically and is shown in Fig. [2j There are two 
parameter regions where the dependence can be essentially simplified, as follows: 

i) for relatively weak field, B ~ 0.1-B e ~ 4 x 10 12 G, the neutrino mean free path can be 
obtained from Eq. fj20[) : 

4.9 m / 219 \ 

A ~ £ .i sine exp Uo.i^i5 sm9 J ' (27) 



6 



E/eV 




10" 2 ID" 1 1 10 1 2 1 3 10 4 B/B e 

Figure 2: The dependence of the neutrino cutoff energy E c on the magnetic field strength for 
the fixed neutrino mean free path value, A = 1 km. 



where B 0A = fl/(0.LB e ), E 15 = E/(W 15 eV), and the cutoff energy corresponding to A = 1 km, 
at -Bq.i = 1) Q = 7r /2, is 



E c ~ 0.4 x 10 17 eV ; (28) 

ii) for relatively strong field, B ~ WB e ~ 4 x 10 14 G, the neutrino mean free path can be 
obtained from Eqs. ([23]), ([SBll : 



3.2 cm / 4.0 \ 

A "ir^ exp v^^p^j' (29) 

where B\q = B/(10B e ), and the cutoff energy corresponding to A = 1 km, at B±q = 1, 9 = tt/2, 
is 

E c ~ 0.6 x 10 15 eV. (30) 

The results obtained show an essential influence of the intense magnetic field on the process 
v — > e~W + width. Despite the exponential character of suppression of the width in a strong 
field, Eqs. (|23|). ([26]) . as well as in a weak field, Eq. (pU]) . the decay width in a strong field is 
greater in orders of magnitude than the one in a weak field, for the same neutrino energy. 



4 Conclusion 

• An influence of a strong external magnetic field on the neutrino self-energy operator is 
investigated. 

• The width of the neutrino decay into the electron and W boson, and the mean free path 
of an ultra-high energy neutrino in a strong magnetic field are calculated. 
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• An energy cutoff for neutrinos propagating in a strong field is defined. The cutoff energy 
at B ~ 5B e ~ 2.2 x 10 14 G is: E c = 10 15 eV . 
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